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Abstract 

We study the ND system by means of a chiral constituent quark model. This model, tuned in the 
description of the baryon and meson spectra as well as the NN interaction, provides parameter-free 
predictions for the charm —1 two-hadron systems. The presence of a heavy antiquark makes the 
interaction rather simple. We have found sharp quark-Pauli effects in some particular channels, 
generating repulsion, due to lacking degrees of freedom to accommodate the light quarks. Our 
results point to the existence of an attractive state, the AD* with (T,S) = (1,5/2), presenting a 
resonance close to threshold. It would show up in the scattering of D mesons on nucleons as a 
D wave state with quantum numbers (T)J P = (l)5/2 _ . This resonance resembles our findings in 
the AA system, that offered a plausible explanation to the cross section of double-pionic fusion 
reactions through the so-called ABC effect. The study of the interaction of D mesons with nucleons 
is a goal of the PANDA Collaboration at the European facility FAIR and, thus, the present 
theoretical study is a challenge to be tested at the future experiments. 

PACS numbers: 14.40.Lb,12.39.Pn,12.40.-y 
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I. INTRODUCTION 



The study of the interaction between charmed mesons and nucleons has become an in- 
teresting subject in several contexts [!}. It is particularly interesting for the study of chiral 
symmetry restoration in a hot and/or dense medium (2J. It will also help in the understand- 
ing of the suppression of the J/\& production in heavy ion collisions [jjj]. Besides, it may shed 
light on the possible existence of exotic nuclei with heavy flavors [4]. Experimentally, it will 
become possible to analyze the interaction of charmed mesons with nucleons inside nuclear 
matter with the operation of the FAIR facility at the GSI laboratory in Germany There 
are proposals for experiments by the PANDA Collaboration to produce D mesons by anni- 
hilating antiprotons on the deuteron. This could be achieved with an antiproton beam, by 
tuning the antiproton energy to one of the higher-mass charmonium states that decays into 
open charm mesons. The D mesons then have a chance to interact with the nucleons inside 
the target material. Since the D mesons are produced in pairs in the antiproton-nucleon 
annihilation process, the appearance of one of those D mesons can be used as tag to look for 
such reactions. These experimental ideas may become plausible based on recent estimations 
of the cross section for the production of DD pairs in proton-antiproton collisions jH[ . Thus, 
a good knowledge of the interaction of charmed mesons with ordinary hadrons, like nucleons, 
is a prerequisite. 

Before one can infer in a sensitive way changes of the interaction in the medium, a 
reasonable understanding of the interaction in free space is required. However, here one has 
to manage with an important difficulty, namely the complete lack of experimental data at 
low energies for the free-space interaction. Thus, the generalization of models describing 
the two-hadron interaction in the light flavor sector could offer insight about the unknown 
interaction of hadrons with heavy flavors. This is the main purpose of this work, to make 
use of a chiral constituent quark model describing the NN interaction |6j] as well as the 
meson spectrum in all flavor sectors 0] to obtain parameter-free predictions that may be 
testable in the future experiments of the PANDA Collaboration. Such a project was already 
undertaken for the interaction between charmed mesons with reasonable predictions |8j, 
what encourages us in the present challenge. 

The paper is organized as follows. In Sec. [Ill we present a description of the quark- 
model wave function for the baryon-meson system, centering our attention in its short- 
range behavior looking for quark-Pauli effects. In Sec. II III we briefly revise the interacting 
potential. Section [IV] deals with the solution of the two-body problem by means of the 
Fredholm determinant. In Sec. [V] we present our results. We will discuss in detail the 
baryon-meson interactions emphasizing those aspects that may produce different results 
from purely hadronic theories. We will analyze the character of the interaction in the 
different isospin-spin channels, looking for the attractive ones that may lodge resonances 
to be measured at PANDA. We will also compare with existing results in the literature. 
Finally, in Sec. I VI I we summarize our main conclusions. 

II. THE BARYON-MESON WAVE FUNCTION 

In order to describe the baryon-meson system we shall use a constituent quark cluster 
model, i.e., hadrons are described as clusters of quarks and antiquarks. Assuming a two- 
center shell model the wave function of an arbitrary baryon-meson system, a baryon Bi and 
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a meson Mj, can be written as: 
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where A is the antisymmetrization operator accounting for the possible existence of identical 
quarks inside the hadrons. In the case we are interested in, baryon-meson systems made of 
iV or A baryons and D or D* mesons, the antisymmetrization operator is given by 



1=1 



(2) 



where Plf T exchanges a pair of identical quarks i and j, and j stands for the light quark 
of the charmed meson. If we assume gaussian Os wave functions for the quarks inside the 
hadrons, the normalization of the baryon-meson wave function ^bm(R) of Eq. ([T]) can be 
expressed as, 



K'mM = K(R) - C(S,T)Af&R) 



(3) 



J\ffa(R) and Af^(R) stand for the direct and exchange radial normalizations, respectively, 
whose explicit expressions are 
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where, for the sake of generality, we have assumed different gaussian parameters for the 
wave function of the light quarks (6) and the heavy quark (6 C ). In the limit where the two 
hadrons overlap (R — > 0), the Pauli principle may impose antisymmetry requirements not 
present in a hadronic description. Such effects, if any, will be prominent for L = 0. Using 
the asymptotic form of the Bessel functions, iL+1/2, we obtain, 
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Finally, the S wave normalization kernel, Eq. ([3]), can be written in the overlapping region 

as 
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where 7 = 1 + -rf-. Thus, the closer the value of C(S,T) to 1 the larger the suppression 
of the normalization of the wave function at short distances, generating Pauli repulsion. In 
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FIG. 1: Norm kernel denned in Eq. ([3]) for L = and four different channels. 



particular, if C(S, T) = 1 the norm goes to zero for R — )■ 0, what is called Pauli blocking [9[. 
C(S,T) is a coefficient depending on the total spin (S) and isospin (T) of the -BjMj two- 
hadron system and is given by 




where the subindices k and i in the exchange operator P ST denote a quark of the baryon 
Bi and a quark of the meson Mj, respectively. \% (Vi) stands for the coupled spin (isospin) 
of two quarks inside the baryon, (Si, T{) are the spin and isospin of the baryon N or A, and 
5*2 is the spin of the meson D or D*. 



TABLE I: C(S,T) spin-isospin coefficients denned in Eq. ((7|). 
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The values of C(S,T) are given in Table [B Similarly to Pauli blocked channels, corre- 
sponding to C(S,T)=1, we will call Pauli suppressed channels those where C(S,T) is close 
to one. We can see that there is one system showing Pauli blocking, this is the AD* with 
(T, S) = (2, 5/2). This is easily understood due to lacking degrees of freedom to accommo- 
date the light quarks present on this configuration. The interaction in this channel will be 
strongly repulsive [9], as we will discuss in Sec. [V] The channel ND* with (T, J) = (0, 1/2) 
where C(S,T) = 2/3 is Pauli suppressed, the norm kernel gets rather small at short dis- 
tances giving rise to Pauli repulsion at short distances as we will see in Sec. [V] We show 
in Fig. [I] the normalization kernel given by Eq. (j3J) for L = and four different channels: 
AD* with (T, J) = (2,5/2), ND* with (T, J) = (0,1/2), ND with (T, J) = (1,1/2), and 
ND* with (T, J) = (0, 3/2). In the first three cases C(S, T) is positive, being exactly one in 
the first case and becoming smaller in the others, what makes the norm kernel to augment. 
In the last case C(S,T) is negative, showing a large norm kernel at short distances and 
therefore one does not expect any Pauli effect at all. 



III. THE TWO-BODY INTERACTIONS 

The two-body interactions involved in the study of the baryon-meson system are obtained 
from the chiral constituent quark model [6|. This model was proposed in the early 90's in 
an attempt to obtain a simultaneous description of the nucleon-nucleon interaction and the 
baryon spectra. It was later on generalized to all flavor sectors |7j. In this model baryons 
are described as clusters of three interacting massive (constituent) quarks, the mass coming 
from the spontaneous breaking of the original SU(2)l <8> SU(2)r chiral symmetry of the 
QCD Lagrangian. QCD perturbative effects are taken into account through the one-gluon- 
exchange (OGE) potential [10]. It reads, 

VOGETO - 4 \ ■ 4 y 2m ^ + 2m 2 + 3m . m J r 2 r .. 4m 2 r 3 I » ' 

where A c are the SU(3) color matrices, = fo//i is a flavor- dependent regularization scaling 
with the reduced mass of the interacting pair, and a s is the scale-dependent strong coupling 
constant given by 0], 

" s(/i) = ^/+\l)hlV (9) 

where a = 2.118, /x = 36.976 MeV and 7 = 0.113 fm -1 . This equation gives rise to 
a s ~ 0.54 for the light-quark sector and a s ~ 0.43 for uc pairs. 

Non-perturbative effects are due to the spontaneous breaking of the original chiral sym- 
metry at some momentum scale. In this domain of momenta, light quarks interact through 
Goldstone boson exchange potentials, 

V x (fij) = Vossifij) + Vbp B (ry) , (10) 
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TABLE II: Quark model parameters. 



m u , d (MeV) 313 
m c (MeV) 1752 
6(fin) 0.518 
6 c (fin) 0.6 
f (MeV fm) 28.170 
a^MeVfm" 1 ) 109.7 
0.54 
3.42 
0.70 
4.2 
230 
0.70 
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g^ h /47r is the chiral coupling constant, Y(x) is the standard Yukawa function defined by 
Y(x) = e~ x /x, = 3 ((7j • fij){fij ■ fij) — <?i • c?j is the quark tensor operator, and H{x) = 
(l + 3/x + 3/x 2 )F(x). 

Finally, any model imitating QCD should incorporate confinement. Being a basic term 
from the spectroscopic point of view it is negligible for the hadron-hadron interaction. Lat- 
tice calculations suggest a screening effect on the potential when increasing the interquark 
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Once perturbative (one-gluon exchange) and nonperturbative (confinement and chiral sym- 
metry breaking) aspects of QCD have been considered, one ends up with a quark-quark 
interaction of the form 



[qflj = nn) => VfcoN^i) + KxjE(r%) + V^(fy) 
[<Mj = era] Vcon^j) + K)ge(^j) 



(13) 



where n stands for the light quarks u and d. Notice that for the particular case of heavy 
quarks (c or b) chiral symmetry is explicitly broken and therefore boson exchanges do not 
contribute. For the sake of completeness we compile the parameters of the model in Table 
ITT1 The model guarantees a nice description of the baryon (N and A) [l2[ and the meson (D 
and D*) spectra Q. Let us also note that the parameters of the model have been tuned in 
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FIG. 2: Different diagrams contributing to the baryon-meson interaction. The vertical solid lines 
represent a light quark, u or d, while the wavy line represents the charm antiquark. The horizontal 
line denotes the interacting quarks. The number between square brackets stands for the number 
of diagrams topologically equivalent. 



the meson and baryon spectra and the NN interaction, and therefore the present calculation 
is free of parameters. 

In order to derive the local B n M m — > B^M\ interaction from the basic qq interaction 
denned above, we use a Born-Oppenheimer approximation. Explicitly, the potential is cal- 
culated as follows, 



V) 



B n M m (L 



ST)^B k M t (L' S'T)(R) - £lS S T T ( R ) ~ £l,ST (°° 



-V S'T 1 



where 



tL 1 s' T 



(R) 



V S'T 
B k Mi 



( R ) I Ei<j=l Vqw ( 



LST 



(R) 



* l b?L( r )\* l b s JJR) 



(14) 



(15) 



In the last expression the quark coordinates are integrated out keeping R fixed, the result- 
ing interaction being a function of the baryon-meson relative distance. The wave function 
m (R) f° r ^ ne baryon-meson system has been discussed in detail in Sec. [Ill We show 
in Fig. [2] the different diagrams contributing to the baryon-meson interaction. The contri- 
butions on the first line are pure hadronic interactions, while the diagrams in the second 
and third lines contain quark-exchanges and therefore will not be present in a hadronic 
description. In Sec. |V] we will analyze in detail the different diagrams to make clear the 
contributions arising from pure quark effects. 



IV. INTEGRAL EQUATIONS FOR THE TWO-BODY SYSTEMS 



To study the possible existence of exotic states made of a light baryon, N or A, and a 
charmed meson, D or D*, we have solved the Lippmann-Schwinger equation for negative 
energies using the Fredholm determinant. This method permitted us to obtain robust pre- 
dictions even for zero-energy bound states, and gave information about attractive channels 
that may lodge a resonance [8J. We consider a baryon-meson system BiMj (Bi = N or A and 
Mj = D or D*) in a relative S state interacting through a potential V that contains a tensor 
force. Then, in general, there is a coupling to the BiMj D wave. Moreover, the baryon- 
meson system can couple to other baryon-meson states. We show in Table III II the coupled 
channels in the isospin-spin (T, J) basis. Thus, if we denote the different baryon-meson 
systems as channel Ai, the Lippmann-Schwinger equation for the baryon-meson scattering 
becomes 

POD 
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x G^E-p^t^f'ip^pp-E) , a,P = A x ,A 2 ,--- , (16) 

where t is the two-body scattering amplitude, T, J, and E are the isospin, total angular 
momentum and energy of the system, £ a s a , ^ 7 s 7 , and IpSp are the initial, intermediate, and 
final orbital angular momentum and spin, respectively, and p 1 is the relative momentum of 
the two-body system 7. The propagators G^(E;p 7 ) are given by 

G,(E ]Pl )= 2/i ; , (17) 

with 

k 2 

E=^, (18) 

where /x 7 is the reduced mass of the two-body system 7. For bound-state problems E < 
so that the singularity of the propagator is never touched and we can forget the it in the 
denominator. If we make the change of variables 

V, = (19) 



where d is a scale parameter, and the same for p a and pp, we can write Eq. ([IE]) as 
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TABLE III: Interacting baryon-meson channels in the isospin-spin (T, J) basis. 
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We solve this equation by replacing the integral from —1 to 1 by a Gauss- Legendre quadra- 
ture which results in the set of linear equations 
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and where w m and x m are the weights and abscissas of the Gauss- Legendre quadrature while 
p 7m is obtained by putting x 7 = x m in Eq. (lT9~j) . If a bound state exists at an energy E B , 

the determinant of the matrix M "^y "'^ (Eb) vanishes, i.e., \M ar ^ TJ (E B )\ = 0. We took 
the scale parameter d of Eq. ([HI as d = 3 fm _1 and used a Gauss- Legendre quadrature with 
TV = 20 points. 



V. RESULTS AND DISCUSSION 

In Figs. [3] and H] we show some representative potentials of the baryon-meson system 
under study. In Fig. [3] we have depicted for two channels the contribution of the several 
pieces of the interacting potential: OGE, OSE and OPE. In Fig. @] we have separated for two 
different channels the contribution of the different diagrams of Fig. [2] in two groups: direct 
terms, those shown in the first line that do not contain quark exchanges, and exchange 
terms, diagrams in the second and third lines of Fig. [2] where quark-exchange contributions 
appear. 

Regarding Fig. [31 there are general trends that have already been noticed for other two- 
hadron systems and we can briefly summarize. For very-long distances (R > 4 fm) the 
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FIG. 3: Contributions to the AD* interaction from the different pieces of the potential in two 
uncoupled channels: (T,J) = (2,1/2) (left panel) and (T,J) = (1,5/2) (right panel). 
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FIG. 4: Direct and exchange contributions to the ND (T,J) = (0,1/2) (left panel) and AD* 
(T, J) = (1,5/2) (right panel) potentials. The dashed-dotted line represents the contribution of 
the direct terms, the dashed line stands for the exchange terms, and the solid line indicates the 
total potential. 

interaction comes determined by the OPE potential, since it corresponds to the longest- 
range piece. The OPE is also responsible altogether with the OSE for the long-range part 
behavior (1.5 fm < R < 4 fm), due to the combined effect of shorter range and a bigger 
strength for the OSE as compared to the OPE. The OSE gives the dominant contribution 
in the intermediate range (0.8 fm < R < 1.5 fm), determining the attractive character of 
the potential in this region. The short-range (R < 0.8 fm) potential is either repulsive or 
attractive depending on the balance between the OGE and OPE. In general, one can say 
that the OGE is mainly repulsive while the OSE is mostly attractive, see the two examples 
in Fig. [31 Thus, in most cases the character of the OPE at short range determines the final 
character of the total potential. 

It is also interesting to analyze the interaction in terms of the different diagrams plotted 
in Fig. [21 as it has been done in Fig. |H The dynamical effect of quark antisymmetrization 
can be estimated by comparing the total potential with the one arising from the diagrams in 
the first line of Fig. [2J which are the only ones that do not include quark exchanges. Let us 
note however that Pauli correlations are still present through the norm, in the denominator 
of Eq. (fl~5l) . To eliminate the whole effect of quark antisymmetrization one should eliminate 
quark-Pauli correlations from the norm as well. By proceeding in this way one gets a genuine 
baryonic potential, that we call direct potential. The comparison of the total and direct 
potentials reflects the quark antisymmetrization effect beyond the one-hadron structure. 
In Fig. H] we have separated the contribution of the direct and exchange terms for two 
different partial waves. As can be seen the direct potential is always attractive, due both 
to the contribution of the scalar exchange interaction and the absence of the one-gluon 
exchange (contributing only through quark-exchange diagrams). However, the character of 
the exchange part, containing quark-exchange diagrams that would therefore be absent in 
a pure hadronic description, depends on the sign of the color-spin-isospin coefficients. The 
sign of the dominant quark-exchange diagram, V34P34, is crucial for determining whether 
the exchange contributions are attractive or repulsive. Thus, quark-exchange diagrams give 
repulsion for the ND (T, J) = (0, 1/2) channel (see left panel of Fig. HJ), but attraction for 
the AD* (T, J) = (1, 5/2) (see right panel of Fig. H|), determining the final character of the 
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A D* (T, J) = (2, 5/2) 
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FIG. 5: Interacting potential in two different channels showing Pauli effects: ND* (T, J) = (0, 1/2) 
(left panel) and AD* (T, J) = (2,5/2) (right panel). See text for details. 



interacting potential. Thus, dynamical quark-exchange effects play a relevant role in the 
ND interaction with observable consequences as we will discuss below, being responsible for 
the repulsive character of the ND interaction at short distances in some partial waves. 

The potential becomes strongly repulsive in those cases where we have Pauli blocked, 
C(S,T) = 1, or Pauli suppressed channels, C(S,T) close to one. In the left panel of Fig. [5] 
we have drawn the ND* (T, J) = (0,1/2) interaction. As seen in Table H] the value of 
C(S,T) = 2/3 is close to one, suppressing the overlapping of the wave function at short 
distances, see Fig[TJ This gives rise to a strong repulsion at short range. The dashed-dotted 
line in this figure represents the direct potential, the one obtained at hadronic level neither 
considering exchange diagrams in the norm kernel nor in the interacting potential. As can 
be seen, the effect of quark-exchanges is rather important in what we have called Pauli 
suppressed partial waves. Hadronic interactions would therefore not be able to account for 
the consequences of Pauli effects as it was already noticed in a comparative study of the NA 



interaction by means of hadronic or quark-based models 13] . In the right panel of Fig. [5] 
we have drawn the AD* (T, J) = (2,5/2) interaction. In this case C(S,T) is exactly one, 
forbidding the overlapping of the two hadrons for R = 0. All contributions are very strong 
at short distances due to the behavior of the norm kernel and the total interaction becomes 
strongly repulsive. These forbidden states would have to be eliminated by hand in the RGM 



treatment of the two-hadron systems [14| . The existence of such a strong repulsion has also 



been observed in the NA system and may be concluded from the NA phase shift behavior 



derived from the nd elastic scattering data [15 



Using the interactions described above, we have solved the coupled channel problem of 



TABLE IV: Character of the interaction in the different baryon-meson (T, J) channels. 





T = T= 1 


T = 2 


J = 


1/2 Repulsive Repulsive 


Weakly repulsive 


J = 


3/2 Weakly attractive Weakly repulsive 


Attractive 


J = 


5/2 Attractive 


Strongly repulsive 



11 




-40 -20 -40 -20 

E(MeV) E(MeV) 



FIG. 6: (T, J) = (0,3/2) (left panel) and (T, J) = (1,1/2) (right panel) Fredholm determinant. 

the baryon-meson systems made of a baryon, N or A, and a meson, D or D* as explained in 
Sec. HYl The existence of bound states or resonances will generate exotic states with charm 
— 1 that could be identified in future experiments of the PANDA Collaboration at the FAIR 
facility [l[ . In Table IIVI we summarize the character of the interaction in the different (T, J) 
channels. It can be observed how all (T, J) channels containing Pauli blocked or Pauli 
suppressed states are repulsive: (2,5/2), (0, 1/2), and (1, 1/2). Thus, the Pauli principle at 
the level of quarks plays an important role in the dynamics of the ND system. 

In Figs. [6] and |8] we have plotted the Fredholm determinant for some representative 
channels. Fig. [6] shows the weakly attractive (T, J) = (0,3/2) and the weakly repulsive 
(T, J) = (1, 1/2) Fredholm determinants. In Fig. [7]we have plotted all potentials contribut- 
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FIG. 7: Baryon-meson potentials contributing to the (T, J) = (1, 1/2) channel. 
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FIG. 8: (T,J) = (2,3/2) (left panel) and (T, J) = (1,5/2) (right panel) Fredholm determinant. 

ing to the (T, J) = (1, 1/2) channel that, as seen in the right panel of Fig. |6l is weakly 
repulsive. The interaction in the lightest two-hadron systems contributing to this channel, 
ND and ND*, has a repulsive character and there is a weak coupling to the attractive but 
heavier two hadron system AD*, giving an overall repulsive channel. Fig. |H] shows the Fred- 
holm determinant for the two attractive channels of the ND system: (T, J) = (2,3/2) and 
(T, J) = (1,5/2). The (T, J) = (1,5/2) is the most attractive one, showing a bound state 
with a binding energy of 3.87 MeV. This channel corresponds to a unique physical system, 
AD*. This situation is rather similar to the one we found in the A A system [16j . predicting 
an S wave resonance with maximum spin. Experimental evidence for such a resonance was 
already reported in the NN scattering data, where the resonance appeared in the 3 D^ partial 




FIG. 9: Binding energy, in MeV, of the (T,J) = (1,5/2) channel as a function of the harmonic 
oscillator parameter used for the charm quark. 
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FIG. 10: (T, J) = (0, 1/2) Fredholm determinant obtained with the direct potentials. The dashed 
line stands for the single channel problem ND and the solid line to the couple channel problem 
using also the ND* state. 



wave, thus as a D wave in the NN system [17] . This prediction has been recently used as a 
possible explanation of the measured cross section of the double-pionic fusion of nuclear sys- 
tems through the so-called ABC effect [18j|. The formation of an intermediate A A resonance 



with the isospin, spin, parity and mass found in Ref. 16j ((T)J P = (0)3 + and M =2.37 
GeV) allowed to describe the cross section of the double-pionic fusion reaction pn — > cItt 7t . 
In the present case, the bound state in the (T, J) = (1,5/2) AD* channel would appear in 
the scattering of D mesons on nucleons as a D wave resonance, what could in principle be 
measured at PANDA. Such a state will have quantum numbers (T)J P = (l)5/2~. Just to 
make sure that our conclusion does not depend on our choice of the parameter of the charm 
quark wave function, that does not come out from the solution of the meson spectra, we 
have calculated the binding energy as a function of b c . The result appears in Fig. [9j showing 
a bound system for any reasonable choice of such a parameter, giving thus confidence to our 
prediction. 



In Refs. [191 . |20| the ND system has been analyzed by means of a hadronic model us- 
ing Lagrangians satisfying heavy quark symmetry and chiral symmetry. They arrive to the 
conclusion that the (T)J P = (0)l/2~ channel is the most attractive one presenting a bound 
state of around 1.4 MeV. To emphasize the importance of quark-exchange effects we have 
repeated for this channel the calculations explained in Sec. [IV] using only the direct poten- 
tials, those without quark-exchange effects. These potentials, that would correspond to a 
purely baryonic interaction, are drawn as a dashed-dotted line in the left panels of Fig. HI 
the (T, J) = (0, 1/2) ND interaction, and Fig. the (T, J) = (0, 1/2) ND* interaction. In 
both cases the importance of quark-exchange effects can be seen. These interactions would 
be the hadronic potentials of an effective theory without quark degrees of freedom. The 
results obtained neglecting the contribution of quark-exchange effects are shown in Fig. [TU1 
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As we can see in both cases, single channel or coupled-channel calculation, a bound state 



appears, in agreement with the conclusions of Refs. [19|, |20] . Once again, this comparison 
makes evident the great importance that quark-exchange effects may have in the system 
under study and it also represents a sharp example of a system where the quark-exchange 
dynamics may have observable consequences. A future effort in the study of the ND system 
will provide us with evidence to learn about the importance of quark-exchange dynamics. 

Let us finally mention that there are recent estimations in the literature |2l) about the 
ND cross section based on hybrid models considering meson-exchanges supplemented with 
a short-range quark-gluon dynamics. The important conclusion of that work is that the 
predicted ND cross sections are of the same order of magnitude as those for the NK system, 
but with average values of 20 mb, roughly a factor two larger than for the latter system. 
Thus, the study of the ND interaction could be a plausible challenge for the PANDA 
Collaboration. 



VI. SUMMARY 



Summarizing, we have studied the ND system at low energies by means of a chiral 
constituent quark model that describes the baryon and meson spectra as well as the NN 
interaction. Due to the presence of a heavy antiquark the interaction becomes rather simple 
and parameter-free predictions can be obtained for the ND system. We have analyzed 
in detail the interaction in the different isospin-spin channels, emphasizing characteristic 
features consequence of the contribution of quark-exchange dynamics. Quark-Pauli effects 
generate a strong repulsion in some particular channels due to lacking degrees of freedom 
to accommodate the light quarks. Such effects have observable consequences generating 
repulsion in channels that otherwise would be attractive in a hadronic description. We have 
traced back our results to the previous analysis of the NA and AA systems with peculiar 
predictions supported by the experimental data. We have found only one bound state in the 
AD* (T, S) = (1, 5/2) system. This state, (T)J P = (l)5/2~, will show up in the scattering of 
D mesons on nucleons asaf wave resonance. Such a resonance resembles our findings in the 
AA system that offered a plausible explanation to the cross section of double-pionic fusion 
reactions through the so-called ABC effect. The existence of this state is a sharp prediction 
of quark-exchange dynamics because in a hadronic model the attraction appears in different 
channels. Finally, it is important to emphasize that theoretical estimations indicate that the 
ND cross section may be attainable in the future facility FAIR and therefore the predicted 
resonance may be a challenge for the study of the PANDA Collaboration. This objective 
may constitute a helpful tool in discriminating among the different scenarios used to describe 
the dynamics of heavy hadron systems. 
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